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Abstract 

We prove some variants of the exponential formula and apply them to the 
multivariate Tutte polynomials (also known as Potts-model partition functions) 
of graphs. We also prove some further identities for the multivariate Tutte 
polynomial, which generalize an identity for counting connected graphs found 
by Riordan, Nijenhuis, Wilf and Kreweras and in more general form by Leroux 
and Gessel, and an identity for the inversion enumerator of trees found by 
Mallows, Riordan and Kreweras. Finally, we prove a generalization of Mdbius 
inversion on the partition lattice. 
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1 Introduction 

Let {c n }™ =1 be a sequence of coefficients, and let 



n! 

n=l 



be the corresponding exponential generating function, considered as a formal power 
series in the indeterminate x. Now let q be another (commuting) indeterminate 1 , and 
define 



A q (x) = exp[qC(x)} = J^aniq)— , (1.2) 



nl 

n=0 

so that A q (x) = A(x) q [we write A(x) = A\(x)\. It is easy to see that: 

(a) Each coefficient a n (q) is a polynomial in q, of degree at most n; indeed we have 
the explicit formula 2 

<-— u ni,. . . ,n e > 1 J— 1 

£ rij = n 

(b) ao(g) = 1. 

(c) For ti^O, the polynomial a n (q) has zero constant term, so that we can write 
a n (q) = qa n (q) where a n (q) is a polynomial of degree at most n — 1; moreover 
we have a n (0) = c n . 

For some purposes it is convenient to introduce the modified generating function 

A q (x) - 1 exp[gC(a;)] - 1 ^ . N x n 



A q (x) = = ~ 1 = 5>(?A ■ (1.4) 

77;. 

n=l 

The coefficients {a n (q)}^ =0 furthermore satisfy the following identities [23] : 3 

^2 (™J a>k(qi) a n - k (q 2 ) = a n (q± + q 2 ) (1.5) 

5^ (I) fc " fc ^ 1 ) a "- fc ^ 2 ) = n «n(^i + ^) for n > 1 (1.6) 
fc=i ^ ' 



1 Let us stress that this g has nothing to do with the q of (/-series; rather, it is connected with the 
q of the g-state Potts model [44,68,69] and of the multivariate Tutte polynomial [57]. 

2 This formula holds also for n = because the integer has exactly one ordered partition 
(ni, . . . , ng) into zero or more strictly positive parts: namely we have £ — parts, and the multino- 
mial coefficient and the empty product both equal 1, correctly giving ao(q) = 1. 

3 See Gould [23] for fascinating comments on the history of identities equivalent to (1.6) [or to 
(1.10) below], which apparently go back to Euler (1748), Rothe (1793), Hindenburg (1796), von 
Ettingshausen (1826) and Hanstead (1881), among others. 
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Indeed, (1.5) follows from writing the trivial identity 



A qi (x)Ag 2 (x) 



(x) 



(1.7) 



and extracting the coefficient of x n , while (1.6) follows analogously from 



A q2 (x) 



x 



d A qi (x) — 1 



dx 



x 



d A qi+q2 (x) — 1 



dx 



qi + q2 



(1.8) 



[23, eq. (4.3)], which is an easy consequence of (1.2). Specializing (1.6) to qi = 0, 
or alternatively applying x-jj^ directly to (1.2), we obtain the well-known recursion 
relation [67, Theorem 3.10.1] 



a n (q) 



k=i 



n — 1 
k- 1 



c k a n ^ k (q) for n > 1 , 



(1.9) 



which allows the {a n (q)} to be calculated given the {c„}, or vice versa. More generally, 
we can form the linear combination n x (1.5) — (gi + g 2 ) x (1-6) to obtain the recursion 
relation 4 



a«(g 2 ) = Yl 

k=l 
n 

= E 



L!) (?i+?2) - D ?i 



fc=i 



n — 1 
fc - 1 



?2 



n- A 



a fc (?i) a„_ fe (g 2 ) (1.10a) 
a*(5i) a„_ fe (g 2 ) (1.10b) 



for n > 1, which allows the {a^fe)} to be calculated given the {a„(gi)} for any fixed 
qi, and which reduces to (1.9) when q\ = 0. We do not know whether there exist 
useful analogues of (1.5) / (1.6) with higher powers of k on the left-hand side. 

Combinatorially, all these identities are versions of the exponential formula [5,13, 
14,60,64,67], which relates the weights c n of "connected" objects to the weights a n (q) 
of "all" objects, when the weight of an object is defined to be the product of the 
weights of its "connected components" times a factor q for each connected compo- 
nent. 5 This idea can be formalized in the following combinatorial model: Let II n be 
the set of partitions of [n] = {1, . . . ,n} into (zero or more) nonempty blocks; give 



4 See [23, eq. (1.3)] or [59, Exercise 1.39, pp. 50 and 62] for the case q\ = 1, which is of course 
trivially equivalent to the general case of qi ^ 0. 

5 In the theory of combinatorial species [5], the exponential formula is merely a special case of 
the more general concept of composition (or substitution) of species F o G [5, section 1.4], which is 
defined so that the exponential generating function of F o G is the composition of the exponential 
generating functions of F and G. When F is taken to be the species E of sets (ensembles in French) 
and H is some species for which there exists a reasonable notion of connectedness, then, letting H c 
be the subspecies of connected -ff-structures, we have the isomorphism of species H = E o H c ; since 
the exponential generating function of E is exp, this immediately implies the exponential formula 
H(x) = cxp(H c (x)). This result also generalizes to weighted and multi-sort species. See [34] for 
results closely related to those of the present paper, formulated in the language of species. 
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t 

a partition n = {ni, . . . , TT e } a weight w(ir) = q E Yl c k;|> where {cfc}^ are arbitrary 

i=i 

coefficients 6 ; and let a n (q) = ^2 w(ir). Then (1.3) is immediate: the multinomial 

7ren n 

coefficient counts the ordered partitions of [n] into blocks of sizes ni, . . . ,ri£, and the 
l/£\ converts this to unordered partitions. As for (1.5)/(1.6), we first observe that 
a n (Qi + 92) can be computed by allowing each block to be colored either "1" (with 
weight q\) or "2" (with weight 52). Identity (1.5) follows by defining S C [n] to be 
the set of elements colored "1", and setting k = \S\. Analogously, a n (qi + q 2 ) can 
be computed by coloring and weighting blocks as before, except that the block con- 
taining one specified element % G [n] (say, i — 1) must be colored "1" and it receives 
weight 1 instead of weight q\. Defining again S C [n] to be the set of elements colored 
"1", and setting k = \S\, we obtain (1.6) in the form 7 

a„(q 1 + q 2 ) = j>] ^ ^ flfc(gi) Qn-fcfe) for ra > 1 . (1.11) 

Sequences of polynomials {a n (q}} satisfying (1.2) [or equivalently satisfying (1.5) 
and not identically zero] have been termed sequences of binomial type by Rota and 
collaborators [12,16,39,48,49,51] and studied by means of the umbral calculus [18, 
49, 52]. 8 The corresponding sequences {a n (q)/n\} have been termed sequences of 
binomial type by Labelle [32] and convolution families by Knuth [30] and Zeng [70]; 
these authors used elementary formal-power-series methods closely resembling those 
used here. A purely combinatorial approach to sequences of binomial type, employing 
the theory of species, has been developed by Labelle [33] (see also [5, section 3.1]). 

Further identities, generalizing (1.5) / (1.6), can be obtained by a powerful trans- 
formation suggested by Knuth [30]. Start with any formal power series A(x) with 
constant term 1; set A q {x) = A(x) g and define polynomials {a n (q)} by (1.2). Now let 
t be a parameter, and define a formal power series A(x; i) by the implicit equation 

A(x;t) = A(xA(x;tf) . (1.12) 

Since A(x; 0) = A(x), we can view the family {A(x; £)}t e R as a special class of "per- 
turbations" of A(x). Now define polynomials {a n (q;t)} by the obvious formula 

00 n 

A(x;ty = ]Ta n (g;t)^-. (1.13) 

r. Til 



6 The coefficients Ck may in turn be generating functions for some additional structure on the set 
[k] — for instance, a graph, a digraph, a hypcrgraph, etc. — under the condition that this structure 
is "connected". 

7 See also [45, Theorem 2.9]; and see [5, p. 188, Exercise 3.1.20(a)] for the special case q\ = 1. 

8 It is not hard to see by induction that if {a n (q)} satisfies (1.5) [or even the special case of 
(1.5) with qi = qz] and is not identically zero, then dega n < n with ao = 1 and a„(0) = for 
n > 1; moreover, with a little extra work it can be shown that A q (x) is necessarily of the form (1.2): 
sec [16, Theorem 4.1] [30,70] [60, Exercise 5.37, pp. 87-88 and 131-132]. 

Let us also remark that Rota et al. [39,51] and many subsequent authors [12,16,48,49], in defining 
"sequence of binomial type", impose the additional condition that dega„ = n exactly (i.e. ci 7^ 0); 
but since this condition is irrelevant for our purposes, we prefer not to impose it. 
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For each t, these polynomials form a family of the type (1.2) and thus satisfy the 
identities (1.5) / (1.6). On the other hand, a straightforward calculation using the 
Lagrange inversion formula [60] yields the remarkable relation 9 

q 

a n (q;t) = — ■ — -a n (q + nt) = qa n (q + nt), (1.14) 
q + nt 

or even more simply 

a n (q;t) = a n (q + nt). (1.15) 

[In particular, a n (q; t) and a n (q; t) are polynomials jointly in q and t, so we can treat 
t as an indeterminate if we wish.] It follows that, for any family {a n (q)} of type (1.2), 
we have the following identities generalizing (1.5) /(1.6): 10 

^2 ( ? j giQfe(gi + kt) q 2 a n -k{q2 + (n- k)t) = (qi + q 2 )a n (qi + q 2 + nt) (1.16) 

^ \V\ ka k (qi + kt) q 2 a n -k(q 2 + (n - k)t) = na n (q 1 + q 2 + nt) for n > 1 
k=i ^ / (1.17) 



where we understand that qao{q) = 1. Multiplying (1.17) by t and summing it with 
(1.16), we obtain the alternate form 

^2 \l) ak(qi + kt) q 2 a n - k (q2 + (n - k)t) = a n (q x + q 2 + nt) . (1.18) 
fc=0 ^ ' 

If in (1.18) we interchange k <-> n — A; and gi <-> g 2 and then replace g 2 by q 2 — nt, we 
obtain the slightly simpler form 

^-■J (fc) + kt ) a n-k{q2 - kt) = a n (qi + q 2 ). (1.19) 

A combinatorial proof and interpretation of (1.16) / (1.18) in terms of the theory of 
species has been given by Labelle [33]. 11 We call (1.16)— (1.19) Abel-type extensions 
of (1.5)/ (1.6) because in the simplest case A(x) = e x , a n (q) = q n they reduce to 
Abel's celebrated extensions of the binomial theorem (also discovered simultaneously 
by Cauchy) [28,46,61]. 

These Abel-type formulae can alternatively be derived as immediate consequences 
of the Pfaff-Cauchy derivative identities (generalizations of Leibniz's rule for the nth 



9 See [60, Exercises 5.37(e) and 5.58, pp. 87-88, 99, 131-133 and 148], where the argument is 
attributed to Eric Rains and Ira Gessel. Please note also that if we write f(x) = x/A{x) t and 
g(x) = xA(x; t)*, then / and g are compositional inverses (as observed in [30,32] for t = 1). 

10 Already Rota et al. [51, p. 711, Proposition 4] and Reiner [45, Theorem 2.10] noticed that, for 
any t, the family {a n (q; t)} defined by (1.14) is of type (1.2). But they do not seem to have noticed 
the relation (1.12), which to our knowledge goes back to Knuth [30]. 

11 See also [5, p. 189, Exercise 3.1.22] for a brief summary. 
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A(x) 


a n (q) 


(1.5)/(1.6) 


(1.16)-(1.19) 


e x 

l+x 
(1-x)- 1 
(l + axf 


q n 

q n N 
q n 1 

a n (/3q)* J 


binomial 
theorem 

Chu-Vandermonde 
convolution 


Abel-Cauchy 
binomial 
identities 

Rothe's generalized 
Chu-Vandermonde 
convolution 


exp(e x — 1) 


n 

n n\q) — u; q 

k=0 

(Bell polynomial [11]) 






exp[— xj (1 — x)\ 


L ( n 1] {q) = EflC:D(-?) fe 

k=0 

(Laguerre polynomial [1]) 







Table 1: Some classical identities that are special cases of (1.5) / (1.6) and (1.16)- 
(1.19). See [21-24,28,46,61] for more information concerning Abel-Cauchy and Chu- 
Vandermonde-Rothe identities. The identity (1.5) is well known for the Bell [11, 
p. 136, eq. 3n] and Laguerre [2, p. 192, eq. 41] [36, p. 96, Problem 20] polynomials. 
We do not know any references for (1.6) or (1.16)— (1.19) in these cases. 



derivative of the product of two functions) recently presented by Johnson [28]. Indeed, 
setting v = A(x) qi , w = A(x) q ' 2 , <fi = A{x) t in [28, eq. (1.1)] yields (1.19), while making 
the same substitution in [28, eq. (1.4)] yields (1.16). 

Some classical examples of families {a n (q)} and the corresponding identities (1.5)/ 
(1.6) and (1.16)-(1.19) are listed in Table 1. 

In this note we begin by proving some easy generalizations of the foregoing for- 
mulae to the case in which the integer index n is replaced by a multi-index n over 
a finite ground set V. We then show some applications of the latter formulae to the 
multivariate Tutte polynomial [57] of a graph G = (V, E). Finally, we derive some fur- 
ther identities for the multivariate Tutte polynomial, which generalize an identity for 
counting connected graphs found by Riordan, Nijenhuis, Wilf and Kreweras [31,40,47] 
and in more general form by Leroux [37] and Gessel [17] and an identity for the in- 
version enumerator of trees found by Mallows, Riordan and Kreweras [31,38]. In 
the appendix we prove an apparently new generalization of Mobius inversion on the 
partition lattice. 

We stress that many (though not all) of the formulae derived in this paper are 
already known. But we think that our approach provides a unified perspective that 
may be of interest. 
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2 Variants of the exponential formula 



Let V be a finite ground set; we shall use multi-indices n = (ni) ie v £ N y and 
commuting indeterminates x = (xi)i e v- We write |n| = ^2 ieV rii an d xI1 — Yliev x T- 
We use the notation n! = n«ey n j' anc ^ analogously for binomial and multinomial 
coefficients. Finally, we denote by (resp. 1) the multi-index with all entries (resp. 
1); and for any subset W C V, we denote by l\y the multi-index taking the value 1 
on W and on V \ W. 

Let R be a commutative ring containing the rationals as a subring; we shall use 
polynomials and formal power series whose coefficients lie in R. (In applications, R 
will usually be either the rationals, reals or complex numbers, or a ring of polynomials 
or formal power series over one of these fields.) One possible approach is to treat the 
coefficients c = (c n ) neN v\{ | as indeterminates, and to take R to be the polynomial 
ring Q[c]; this approach has the advantage of exhibiting the polynomial dependence 
on the coefficients c n . But we need not commit ourselves to any specific choice of the 
ring R; all of our identities will be valid in complete generality. 

So let {c n } neN v be a sequence of coefficients in the ring R, with c = 0; and let 

n 

be the corresponding exponential generating function, considered as a formal power 
series in the indeterminates x. Now let q be another (commuting) indeterminate, and 
define 

A,(x) = exp[gC(x)] = J> n (g)^. (2.2) 

n 

Using the Taylor series for exp, we deduce immediately that 

(a) Each coefficient a n (q) is a polynomial in q, of degree at most |n|; indeed we 
have the explicit formula 12 

«.<«> = E£ E !„.." ,JII («> 

£=0 m,...,n^O V ^ ' l/ j=l 
E"j=n 

(b) a (q) = 1. 

(c) For n^O, the polynomial a n (q) has zero constant term, so that we can write 
a n(<?) = Qttniq) where a n (q) is a polynomial of degree at most |n| — 1; moreover 
we have a n (0) = c n . 



12 This formula holds also for n = because the multi-index has exactly one ordered partition 
(ni, . . . , n^) into zero or more nonvanishing parts: namely, we have I = parts, and the multinomial 
coefficient and the empty product both equal 1, correctly giving ao(q) = 1. This same observation 
applies to all subsequent formulae of a similar type, such as (2.7) ff. 
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Conversely, it is not hard to see that if (c n ) neN v\{ } is any sequence and we define 
(a n (q)) neN v by (2.3), then the latter family satisfies (2.2). 

Using qC(x) = log A q (x.) together with the Taylor series for log(l + z), we obtain 
the inverse formula 

^ = e^ e ( ni ,. n ,Jn%M («) 

1=1 m,...,n^O V X ' ' lJ j=l 



E 



n,=n 



or equivalently 



C n = 



f i t r~ E ( ni , n , n ,)n%(«)- (-) 

1=1 ni,...,n^O x ' j=l 

E n 3= n 

It is a nontrivial fact that the right-hand side of (2.5) is independent of q. 

Now let r be another indeterminate, and let us use the identity A rq {x) = A q (x) r 
together with the Taylor series 

(l + z) r = exp[rlog(l+E)] =YJi\ Z ' (2 ' 6) 

t=o l - 

where r- — r(r — 1) ■ ■ ■ (r — £ + 1) denotes the falling factorial. (If we like, we can 
alternatively use the notation Q = r-/£\, which is a well-defined polynomial in the 
indeterminate r.) We obtain 

«»w = E^ E ( ni ,. n .,JlK<«>- < 2 - 7 > 

t=0 ni,...,n^O W ' ' £/ j=l 
Enj=n 

In view of the identity (— r)- = (— l) e r e where r 1 = r(r + 1) • • • (r + £ — 1) = (r+€ — 1)- 
denotes the rising factorial, we also have 

-h^E^ E ( ni , n , n ,)n^w- <«) 

£=0 ni,...,n^0 V X ' ' £/ j=l 

E n 3= n 

Since both sides of the identities (2.7)/(2.8) are polynomials in r, we can specialize 
these identities to r integer or rational, or more generally to r being any element of 
the ring R. If r is a positive integer, the sum in (2.7) can obviously be restricted to 
£ < r. Slightly (but not much) less trivially, if r is a positive integer, then by using 

t (r + l-1)! (e+lf* 

£\ £\{r-l)\ (r-1)! 1 ' ' 
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we can rewrite (2.8) as 

<?=0 V ^ ni,...,n^O W ' ' l/ j=l 

in which the coefficient (£ + l) r_1 /( r ~~ 1)' i s a polynomial of degree r — 1 in £. In 
particular, for r = 1 we have 

«.(-«) = E(-D f e L " J n «■*<»>■ f 2 - 11 ) 

^=0 ni,...,n^0 V X ' " ' ' £/ J=l 

E"3=n 

By formally setting r = q^Azi (this is justifiable if we work in a ring of formal 
Laurent series in the indeterminate q 1 ), we obtain the identity 

= Ejrfn^-^)) e L n Jn^(9i). ( 2 - 12 ) 

E n j= n 

which is nice because both sides are manifestly polynomials in the indeterminates q\ 
and q 2 . Alternatively, for n ^ we can pull out the j — factor and write 

«•(») = E^(n(»-i«)) E L " J ■ (2-13) 

Enj=n 

Conversely, it is not hard to see that if qi is any element of the ring R and (a n (9i)) n eN v ' 
is any sequence, and we define (a n (q)) n€N v by (2.12), then the latter family satisfies 
(2.2). 

The following special cases of (2.12)/(2.13) are of particular interest: 

(a) For qi = 0, the identity (2.12) reduces to (2.3). 

(b) For q 2 = 0, we obtain the inverse formula (2.4)/(2.5). 

(c) For q 2 = —qi, we obtain the closely related formula (2.11). 

It is easy to show that the coefficients {a n (<?)}neN v satisfy the identities 

Qk(gi) Qn-k(g2) = a n (q 1 + q 2 ) ( 2 - 14 ) 
E (k) Q n-k(?2) = nia n (qi + q 2 ) for n 7^ (2.15) 
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Indeed, (2.14) follows by expanding out the identity A qi (x)A q2 ('x) = A qi+q2 (x) and 
extracting the coefficient of x n , while (2.15) follows analogously from 



A q2 (x) 



' d AJx) - 1 

Xi r% 

oxi qi 



d A qi+q2 (x) - 1 
dXi qi + q 2 



which is an easy consequence of (2.2). In particular, specializing (2.15) to q± = 0, we 
obtain the recursion relation 

a n(q) = q ^ C?_ *') c k a n -k(g) whenever rii > 1 , (2.17) 

which allows the {a n (q)} to be calculated given the {c n }, or vice versa. [Here 6i is 
the vector with entry 1 at element i and elsewhere.] More generally, we can form 
the linear combination rii x (2.14) — (gi + q 2 ) x (2.15) to obtain the recursion relation 



Ok(gi) a n -k(g , 2) whenever n { > 1 



(2.18) 

which allows the {a n (g 2 )} to be calculated given the {a n (qi)} for any fixed q 1 , and 
which reduces to (2.17) when q\ = 0. 

We do not know whether there exist useful analogues of (2.14)/(2.15) with higher 
powers of k on the left-hand side. 

Finally, we can deduce the Abel-type extensions of (2.14)/(2.15): 

^ ( k ) qi ^ qi + k ' <?2a n -k(<?2 + (n - k) • t) = (q x + q 2 )a ll (q 1 + q 2 + n • t) 
k V J (2.19) 

^ \ k ) ki ^ qi + k ' *) ^a n _ k (g2 + (n - k) • t) = m a n (qi + q 2 + n • t) for n ^ 
k^o V / (2>20) 

where t = (ti)i S y are commuting indeterminates and we understand qa (q) = 1. 
These formulae can be proved by repeated use of the Lagrange inversion formula, 
once in each variable. For details we refer to the paper of Zeng [70, pp. 225-226]. For 
further information on classical specializations of (2.19)/(2.20), see [42,43,58,61,70]. 

Remark. We wonder whether (2.19)/(2.20) might have an extension to identities 
involving a not-necessarily-diagonal matrix T = (tij)ij e v of indeterminates, based on 
considering the implicit equation 

e C(x;T) = A ( e C(x;T)T x ) ? (2 .2l) 

which generalizes (1.12). If true, this could be a powerful extension. See Strehl [61] 
for possibly related work in which a matrix of indeterminates is employed. 
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3 Application to the multivariate Tutte polyno- 
mial 



In this section we apply the formulae of Section 2 to deduce some identities for 
multivariate Tutte polynomials. 



3.1 Definitions and basic properties 

Let G = (V, E) be a finite undirected graph with vertex set V ^ and edge set 
E; loops and multiple edges are allowed unless explicitly stated otherwise. Then the 
multivariate Tutte polynomial [57] of G is, by definition, the polynomial 

Z G (q,v) = J2v k(A) H v ^ 

ACE eeA 

where q and v = (t> e ) ee _E are commuting indeterminates, and k(A) denotes the number 
of connected components in the subgraph (V, A). It is also convenient to pull out one 
factor of q by defining 

Za(q,v) = E^ M II^; (3-2) 

ACE eeA 

this is still a polynomial since k(A) > 1 for all A. Finally, let us make the convention 
that if G = (the graph with empty vertex set and empty edge set), then Z = 1 
and Z is undefined. 

If we specialize to v e — v for all edges e, we obtain a two-variable polynomial 
Z G (q, v) that is essentially equivalent to the classical Tutte polynomial [57, sec- 
tion 2.5]. 

Note that at q = 1 we have the trivial formula 

Z G (l,v) = Z G (l,v) = ]J(l + v e ) . (3.3) 

e£-E 

More interestingly, at q = we have the specialization 

Zq(0,v) = C a (v) = 11^, (3-4) 

ACE eeA 
k{A) = 1 

i.e. the generating polynomial of connected spanning subgraphs of G. 

In statistical physics, Z G (q, v) is known as the partition function of the g-state 
Potts model [44,68,69] in the Fortuin-Kasteleyn representation [15,26,29]. This arises 
by virtue of the following identity [15,29]: 

Theorem 3.1 (Fortuin-Kasteleyn representation of the Potts model [15,29]) 

For integer q > 1, we have 



Z G (q,v) = e n 

a: V->[q] e=ij&E 

where 5 denotes the Kronecker delta. 



l + v e 8(a(i),a(j)) 



(3.5) 
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Proof. On the right-hand side of (3.5), expand out the product over e G E, and let 
A C E be the set of edges for which the term v e 5{cr{i), cr(j)) is taken. Now perform 
the sum over maps a: V — > [q\: in each connected component of the subgraph (V, A) 
the "color" a(i) must be constant, and there are no other constraints. Therefore, the 
right-hand side equals 

ACE eeA 

as was to be proved. □ 

Specializing to v e = —1 for all edges e, we obtain the Birkhoff- Whitney [4,66] 
expansion for the chromatic polynomial of G: 

Corollary 3.2 For integer q > 1, the number of proper q- colorings of G is -Pg(<?) = 
'/-esq. I!- 

3.2 Some generating functions 

For the remainder of this section, let us assume that G is a loopless graph. For 
notational simplicity it is convenient to assume also that G has no multiple edges; 
the trivial changes to allow multiple edges can be left to the reader. 

For any multi-index n G N v , let us define G[n] to be the graph obtained from 
G by expanding each vertex i to an independent set consisting of Hi vertices. That 
is, the vertices of G[n] are pairs (i,a) with i £ V and a G [rij], and the edges of 
G[n] are all pairs ((i,a), (j, (3)) with ij G E, a G [rij] and /3 G [rij]. Given a set of 
weights v = (v e ) ee E associated to the edges of G, we assign the weight % to each 
edge {(i,a), (j,/3)) in G[n}. 

In a similar way, let us define G'[n] to be the graph obtained from G by expanding 
each vertex i to a clique consisting of n« vertices. That is, the vertex set of G'[n] is 
the same as that of G[n], and the edges of G'[n] consist of those of G[n] together 
with all pairs ((i, a), (i, (3)) with % G V, a, (3 G [rij] and a ^ (3. Given a set of weights 
v = (ve)e£E associated to the edges of G and another set of weights w = (wj)j g y 
associated to the vertices of G, we assign the weight Vij to each edge ((i,a), {j,(3)) 
(i 7^ j) in G'[n] and the weight Wi to each edge {(i,cx), (i,/3)) in G'[n]. 

Note in particular that if n = l w for some subset W C V, then G[n] = G'[n] = 
G[W], the induced subgraph of G on W. 

The following "master formula" generates the multivariate Tutte polynomials 
(with weights as assigned above) of all the graphs G[n] and G'[n]: 

Theorem 3.3 Let G = (V, E) be a simple graph, and let v = (v e ) e <z E and w = 
(wi)i e v be indeterminates. We then have the following exponential generating func- 
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tions: 



E z °m(?.v)^ = (e ( IK 1 + ««)"■"') (37) 

n \ n ij€i? / 

E^(*'.»)j = fE (IK 1 + %>"•"') (IK 1 +--)"'""" ,/2 )Sy 

n \ n ijeE i£V J 

(3.8) 

where all sums run over n e 

We remark that the expression in large parentheses on the right-hand side of (3.7) 
[resp. (3.8)] is the grand partition function for a lattice gas on the graph G [resp. on 
the graph G° obtained from G by adjoining a loop at each vertex] in which arbitrary 
nonnegative integer occupation numbers n = {ni) ie y are allowed, with fugacities x on 
the vertices and two-particle Boltzmann weights 1 + % on edges ij [and two-particle 
Boltzmann weights 1 + w, on the loops]. See e.g. [53,54] for definitions concerning 
lattice gases. 

When G = Ki (the graph with one vertex and no edges), (3.8) reduces to the well- 
known [6, 17, 19,63,65] exponential generating function for the Tutte polynomials of 
the complete graphs K n : 13 

oo n / oo n \ 1 

E z .fc4 = Ed+")" ( "- 1)/2 ^ • ' 3 - 9 > 



n=0 \n=0 



When G = K 2 , (3.7) reduces to the exponential generating function for the Tutte 
polynomials of the complete bipartite graphs K nitn2 , 

£ z «M^h = \ £ < 1+ "»"'"' i^h.) • (3 - 10) 

ni,n 2 =0 1 z \ni,ji 2 =0 1 z / 

which in turn specializes for v — — 1 to the well-known [60, Exercise 5.6, pp. 73 
and 107-108] exponential generating function for the chromatic polynomials of the 
complete bipartite graphs 14 , 

m,n 2 =0 1- z ' 

More generally, when G — K r , (3.7) gives the exponential generating function for the 
Tutte polynomials of the complete r-partite graphs K nu _ jTlr . 



13 Unfortunately, this formula is usually written in terms of the classical Tutte polynomial Tq(x, y), 
which is related to Zo(q, v) by a change of variables [57, section 2.5] that obscures the very different 
roles played by q and v. 

14 (3.11) can alternatively be derived from the explicit expression [62] for Pif „ 2 (<?) by using a 
well-known Stirling-number identity [25, eq. (7.49)] together with the binomial formula. 
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Note also that if we only want to know ZQ\ n ](q,v) for n < 1 — that is, we want 
to know Zg[w](q, v) for subsets W C V — then it suffices to get correct the terms 
n < 1 on the right-hand side of (3.7); all the other terms can be dropped or altered 
arbitrarily. It follows that 

^(^^^^(^(na+^r)^) • (3.i2) 

\n<l ijeE ' J 

Note that the expression in large parentheses is now the grand partition function for 
a lattice gas with hard-core self-repulsion, i.e. at most one atom can occupy each site. 
In the special case v = — 1, this becomes 

P G[ w](q) = [x 1 -]/ G (x)^, (3.13) 

where Ig( x ) is the multivariate generating polynomial for independent sets of vertices 
in G. The formula (3.13) was found previously by Lass [35, Proposition 3.1]. 

Since the formula (3.7) for G[n] is just a specialization of the formula (3.8) for 
G"[n], obtained by taking w = 0, it suffices to prove the latter. 

First Proof of Theorem 3.3. Since the coefficient of x n on each side of (3.8) is 
a polynomial in q, it suffices to prove the identity for infinitely many values of q in the 
ring R; in particular, it suffices to prove it for all positive integers q. By Theorem 3.1, 
for integer q > 1 we have 



Z G>] (q,v,w) = II 

a: V(G'[n])-»[g] rseE(G'[n\) 



1 + v rs 5(a(r),a(s)) 



(3.14a) 



e n n 

a: V(G'[n])-M \ijeE a g [„,] 

fi e K-] 

n n 

a ,pe [m] 

a</3 



l + v i:j 5(a(i,a),a(j,(3)) 



1 + Wi5(a(i,a),a(i,(3)) 



x 



(3.14b) 



Now, because of the symmetries in the construction of G"[n], the summand in (3.14b) 
depends only on the numbers {ni tT }i e v,re[q] defined by 



n itT = #{ae[ni]: cr(i,a)=r}, 



(3.15) 



and we have 

^G'[n](g,V,w) 



e n 



i n i,r} 
Er=l = n » V * 



n; 



n no-^)"" 1 "" 7 



X 



n n^+wi) niAni ' T ~ i) ' 2 



(3.16) 



K i & V r=l 
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Now multiply by x n /n! and sum over n; the sums over the rii iT now become unre- 
stricted, and the sums for the different values of r decouple; we end up with the 
product of q terms each of which is 

e (n( i+ %) nin (n( i+ ^) niK ~ i)/2 ) • ( 3 - i? ) 

n ij£E ieV 

□ 

It is instructive to give an alternate proof of Theorem 3.3 in which q is treated 
directly as an indeterminate: 

Second Proof of Theorem 3.3. We begin by proving the special case (3.9). 
From the binomial formula we have 

oo n \1 oo b / oo n \ b 

n=l ' / 6=1 \n=l ' / 

= ! + Ef e fSrK i +")'" , "'- i,/2 - 

6=1 ' m,...,n b >l 11 j=l 

(3.18b) 

Now write n = Yl'j=i n ji an d consider rii, ...,7i& as the block sizes in an ordered 
partition n = (tti, . . . , n b ) of [n] into b nonempty blocks. The number of ordered 

6 

partitions with those block sizes is precisely the multinomial coefficient n\j Yl so 

j'=i 

(3.18) is equal to 



1 + £^£ 



n=l vren; 1 1 3=1 



1+u )|fil(|Ti|-l)/2 ? ( 319 ) 



where 11^ denotes the set of ordered partitions of [n], and | vr | denotes the number of 
blocks in the partition ir. Passing to unordered partitions absorbs the factor l/|7r|! 
(as these are partitions of a labelled set), and we get 



oo 

X 



1 + E —\ E ^II^ 1 + ^) kjl(kjhl)/2 • (3.20) 
n=i ' 7ren n j=i 

The factor (1 + u^j'KK'I -1 )/ 2 corresponds to a sum over (simple undirected) graphs 
on the vertex set ttj, with a weight i> for each edge. We can now reinterpret (3.20) 
as a sum over (simple undirected) graphs H on the vertex set [n], with a weight 
V \ E ( H )\^ together with a sum over partitions n E U n that are compatible with if in 
the sense that each block of n is a union of (vertex sets of) connected components of 
H. If H has r connected components, then these can be grouped into b nonempty 
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blocks in ways (where {£} is a Stirling number of the second kind); and it is well 
known [25, eq. (6.10)] that 

E {[) V-=<l r - (3-21) 



6=1 



Therefore, each graph H gets a weight q k( - H ^v^ E( - H ^, and the sum over graphs gives 
Z n (q,v). This proves (3.9). 

The more general formulae (3.7)/(3.8) are derived by a variant of this proof in 
which sums over integers n G N\{0} are replaced by sums over multi-indices n G 
{0}. In the one-dimensional case, we were considering partitions of [n], a collection 
of n distinct objects of a single type. We are now working with multi-indices in N v , 
and so we will be partitioning collections of objects of \V\ different types. We begin 
with some notation. Given n = (ni) ie v G N y , we write [n] for the set of ordered pairs 
{(i,a): i G V and a G [n«]}. Thus |[n]| = J2iev Ui - We sna U think of the subset 
[n]j := {(i,a): a G [rij]} as a copy of [rij] with type i; thus, [n] is a collection of 
|n| objects, with rii objects of type i for each i G V. Note that G'[n] is simply the 
complete graph on the vertex set [n]. Now, given a partition 7r = ( y 7r 1 , . . . ,ir k ) of 
[n] (with all ttj nonempty), we write |-7r| = k. The block sizes B(iZj) are elements of 
N v \ {0} defined by 

B(tt j -) = (ttj n [nJOiev ■ (3.22) 

Thus B(7Tj) lists the number of objects of each type contained in the j'th block of 7r. 
Now let us define the weight 

W(n) := ( Y[(l + v u ,) nin A m(l + ^) n<(n< " 1)/2 ) • (3-23) 
\u'eE ) Kiev ) 

The right-hand side of (3.8) is then 

1 + 5><»>£Y = '+E^fe»1n)^ (3.24a) 



b\ \^ v n! 

n^O / b>l Vn^O 

q- ^ x ^ nj 
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6>1 ' ni,..., n( ,^0 11 J ' j=l 

Now let n = Y^!j=i n j> an d consider ni, . . . , nj, as the block sizes in an ordered partition 
7T of [n] into b nonempty blocks. The number of ordered partitions with those block 

b 

sizes is the (multidimensional) multinomial coefficient n!/ Y\ n yi so (3.24) is equal to 

3=1 

i+E^ESri^w). < 3 - 25 > 

n^O 7T jr' = l 

where the second sum is over ordered partitions of [n]. Passing to unordered parti- 
tions, this becomes 



i + E^E^ M n^ B ^))- ^ 



X 

n 

n^O 7r j=l 
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The factor W(B(iTj)) corresponds to a sum over graphs with vertex set 7Tj, such that 
edges within [n]j have weight Wi and edges between [n]j and [n]^ have weight vu>. As 
before, we can reinterpret (3.26) as a sum over graphs H with vertex set [n], with 
these same edge weights. The final step of the argument — summing over partitions 
7r that are compatible with H — is identical to that in the one- dimensional case: it 
makes no difference whether the vertex set is [n] or [n]. We therefore see that each 
graph H gets weight q k ( H ^ times the product of edge weights, and we are done. □ 



3.3 Consequences for Zc{q, v) 

In view of Theorem 3.3, we can apply all the results of Section 2 with the identi- 
fications 

a n (q) = Z G , [n] (q,v,w) (3.27a) 
a n (q) = Z G , [n] (q,v,w) (3.27b) 

Cn = C G #[ n ](v,w) (3.27c) 

an(i) = ( n( i+ %r n M (n( i+ ^) n<(n<_i)/2 ) ( 3 - 27d ) 

\ijeE J \ieV J 

We refrain from writing out all the formulae, which are obtained by substituting 
(3.27a-d) into (2.3), (2.4)/(2.5), (2.7)-(2.13), (2.14), (2.15), (2.17), (2.18), (2.19) and 
(2.20). Let us simply show a few important results that arise from the partition 
formula (2.7)/(2.12) and from the convolution formula (2.14)-(2.15). 



3.3.1 Partition formulae 

By specializing (2.7)/(2.12) to n = 1, we can obtain a partition formula that 
expresses the multivariate Tutte polynomial of G at q = q2 in terms of the multivariate 
Tutte polynomials of induced subgraphs of G at q — qi, for arbitrary choices of q± 
and q 2 . 

Proposition 3.4 Let G = (V, E) be a finite loopless graph. Then 

Z G (q 2 ,v) = ^(g 2 /#nV](?i' v ). ( 3 - 28 ) 

7ren(v) sevr 

where the sum runs over (unordered) partitions n of the vertex set V , the product 
runs over blocks B ofn, and G[B] denotes the induced subgraph of G on the vertex 
set B. 15 Equivalently we can write 

/W-i \ 

Z G (q 2 ,v) = [n^-jqi) ) Y[ Z G[m ( qi ,v) (3.29) 

7Tgn(V) V j=Q / B£lT 

15 Readers who dislike dividing by indeterminates can reinterpret this formula by writing r = 52/91 
and hence <?2 = rqi, as in (2.7). 
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or 

/M-i \ 

Z G (q 2 ,v) = [n(Q2-jqi))Y[Z G[Bi ( qi ,v). (3.30) 

ne!l(V) \ j=l J Ben 

Proof. Specializing (2.7) [or equivalently (2.12)] to n = 1, we obtain a sum over 
ordered partitions 7r = (tti, . . . , 7r^) of V into £ nonempty blocks. Passing from ordered 
to unordered partitions, we get a factor £\ that cancels the l/£\ in (2.7)/(2.12). □ 

Proposition 3.4 is very powerful because we are free to choose qi as we please; then 
Zg(Q2,v) can be written in terms of Zg[b](qi, v )- The two most important special 
cases are q 1 = and q 1 = l. 16 

Case <7x = 0. Taking the limit q\ — > in (3.28), or equivalently just setting 
q 1 = in (3.29), we obtain 

Z G (q,y) = E ^ II ^(v). (3.31) 

7ren(y) Ben 

[This formula can alternatively be derived from (2.3) by the same reasoning as used 
in proving Proposition 3.4.] Of course, (3.31) also has an obvious direct combinatorial 
proof, based on the definition (3.1): first we classify subsets ACE according to the 
partition ix of the vertex set V that is induced by the connected components of the 
subgraph (V, A); then we sum over ways of connecting up each component. 

Case q 1 = 1. This special case is slightly less obvious, and it is worth stating it 
explicitly: 

Corollary 3.5 Let G = (V, E) be a finite loopless graph. Then 

7ren(y) Ben eeE(G[B]) 

where the sum runs over (unordered) partitions it of the vertex set V , the outermost 
product runs over blocks B of it, and the innermost product runs over edges in the 
induced subgraph G[B}. In particular, 

Caiy) = E (-l) |7rhl (N-l)! II II ( 3 - 33 ) 

7ren(y) Ben eeE(G[B\) 

zq(-i,v) = e (-i) w H! n n ( 3 - 34 ) 

7ren(v) Ben e eE(G[B\) 

where \n\ denotes the number of blocks in it. 

Specializing (3.32)-(3.34) to v e = —1 for all edges e, we have: 



16 The case q\ = — 1 is also of interest, as it is related to acyclic orientations. See Lass [35] for 
some interesting material that may be relevant in the present context. 
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Corollary 3.6 Let G = (V,E) be a finite loopless graph with \V\ = n. Then: 



n 



Pa(q) 



(3.35) 



k=i 



n 



1) 



E( 



l)"- 1 (k- 1)1 Q G (k) 



(3.36) 



k=i 



n 



Pg{ 



1) 



E( 



l) k k\ Q G (k) 



(3.37) 



k=i 



where Qg(^) denotes the number of partitions of V into k nonempty independent 
subsets. 

Of course, (3.35) is well known and has a trivial direct proof by counting proper 
g-colorings; and (3.36)/(3.37) are immediate corollaries of (3.35). 

Let us conclude by giving two direct combinatorial proofs of the fundamental 
Proposition 3.4: one using the subgraph representation (3.1), and the other using the 
coloring representation (3.5). 

Second Proof of Proposition 3.4. Use the partition formula (3.31) [which 
itself had a simple direct combinatorial proof] to expand Z G {q 2 , v) on the left-hand 
side of (3.28) and to expand each Zg[b](qi, v) on the right-hand side. We can then 
interpret the double sum on the right-hand side as a sum over partitions it of V, with 
weight q\ for each block, together with a sum over ways of grouping these blocks into 
nonempty groups, with a weight (q 2 /qi)- if we have k groups. Using the identity 



where n — \ir\, we get exactly the expression on the left-hand side. □ 

It will not escape the reader's notice that the argument using (3.38) is identical 
to the one given earlier using (3.21). 

Third Proof of Proposition 3.4. Since both sides of (3.28) are polynomials in 
qi and q2, it suffices to prove the identity for all pairs of positive integers qi, q2 such 
that r = q%/qi is also an integer. We shall do this using Theorem 3.1. Consider a 
(not-necessarily-proper) coloring a: V — > [q 2 ]. Fix a bijection [q 2 ] — [qi] x [r] and 
write a = (r,ip) where r: V — > [qi] and tp: V — > [r]. We then trivially have 




(3.38) 




a: V^[q 2 ] e=ijeE 




(3.39a) 



E E II i + v e S(r(i),T(j))5(i;(i)^(j)) . (3.39b) 



V>: V-»[r] r: V^[ qi ] e=ij€E L 
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The map ip induces a partition tt^, of V into its nonempty color classes, and the 
remaining sum over r gives (by Theorem 3.1 again) n^e^ ^g[b](qi, v )- On the other 

hand, each partition tx of V arises in this way from different colorings tfj, since 
we have r choices to color the lexicographically first block of tc, r — 1 choices for the 
lexicographically second block, and so forth. This proves (3.28). □ 

Important note: The formulae in this subsection — notably (3.33) together with 
the q — 1 case of (3.31) — are strongly reminiscent of Mobius inversion on the lattice 
of partitions II(V). This is not an accident. In fact, the formula (3.30) is a special 
case of a u qi~q 2 generalization" of Mobius inversion that we present in Appendix A. 
We are grateful to an anonymous referee for drawing our attention to the connection 
of these formulae with Mobius inversion and for challenging us to find the underlying 
general principle. 



3.3.2 Convolution formulae 

In a similar way we can obtain convolution formulae for the multivariate Tutte 
polynomial by specializing (2.14)-(2.15): 

Proposition 3.7 Let G = (V, E) be a finite loopless graph. Then 

Z G {qi + q2,v) = Z G[w](qi,v) Z G[v \ w] (q 2 ,v) . (3.40) 
wcv 

Moreover, for each i G V we have 

Z G (q 1 + q 2 ,v) = Y Z G[W](qi, v) Z G[v \ w] (q 2 , v) (3.41) 
wcv 

W 9 i 



and in particular 



Z G (q,v) = Y qC G[W ](v) Z G[v \ w] (q,v) . 



(3.42) 



wcv 

W 9 i 



We also have 



\V\Z G (q 2 ,v) = [(n + teW\- H\V\ 

es+wcv 



Z G [w](qi: v) Z G [ V \ W ](q 2 , v) . 



(3.43) 



Proof. (3.40) and (3.41) are just (2.14) and (2.15), respectively, specialized to n = 1. 
Further specializing the latter to q± = 0, we obtain (3.42). Finally, (3.43) is obtained 
by subtracting (q 1 + q 2 ) x (3.41) from (3.40) and summing over i G V . □ 

For completeness, let us give a simple direct proof of (3.41). [(3.40) is similar but 
easier, and is left as an exercise for the reader.] 



20 



Alternate Proof of Proposition 3.7. From the definition (3.2) we have 



Z G ( qi + q 2 ,v) = J2(Qi + Q2) k{A) - 1 l[v e . 



(3.44) 



ACE eeA 



The factor (qi + q2) k( - A ^ 1 can be handled by choosing, for each connected component 
of (V, A) other than the component containing the distinguished vertex i, to color it 
either "1" or "2", with a corresponding factor qi or q 2 , and summing over all such 
choices. Now define W C V to be the vertex set corresponding to the union of the 
components colored "1" together with the component containing i. By construction, 
W is compatible with A in the sense that A has no edges connecting W to V \ W . 
Moreover, every W 3 i compatible with A is obtained exactly once by some choice of 
colors "1" and "2". On the other hand, the right-hand side of (3.41) is given precisely 
by such a sum over compatible pairs (W, A) [the ~ on Zg[w] (Qi , v ) accounts for the 
fact that the component containing i carries no color and hence no factor □ 

We do not know whether the Abel-type extensions (2.19)/(2.20), when applied to 
the multivariate Tutte polynomial, are of any interest. 

3.4 Specialization to the complete graphs 

Let us now take G to be the complete graph K n , with all the edge weights v e set 
equal to the same value v. We write Z n (q, v), Z n (q, v) and C n {v) for the corresponding 
polynomials. 

3.4.1 Partition formulae 

Specializing (3.28), we obtain the general identity 



where the sum runs over (unordered) partitions it of [n], say n = {ni, . . . ,ne}, and 
| VI- j denotes the number of blocks in ir (i.e., | vr j = I). Special cases of (3.45) arise at 
particular values of qi. 

Case gi = 0. Specializing (3.31), we obtain the trivial identity 




(3.45) 




(3.46) 
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Case gi = 1. Specializing (3.32)-(3.34), we obtain the less trivial relations 

Z n (q,v) = J2 <l-( l + v ) n(n ~ 1)/2 ~ M (3-47) 



C n {v) = ^(-l)l^~ 1 (|7r|-l)!(l+t;) ri ("- 1 )/ 2 -^ll (3.48) 
7ren„ 

Z n (-l,v) = ^(-l) 1 " 1 |7T|! (l +v )n(n-l)/2~ M (3^) 



TreUr, 



where 



7T = 



E = \ \ n2 - Ew 2 ) ( 3 - 5 °) 



2 , 

l<i<j<£ \ i=l 



denotes the number of "cross-edges" in 7r, i.e. the number of edges in the complete 
^-partite graph with vertex classes 7Ti, . . . , iTg. Equivalently, 

n{n ~ l) - h\\ = E 1 *" 11 ( 3 - 51 ) 

i=i 

is the number of "internal edges" in ir. The formulae (3.47)-(3.49) will play an im- 
portant role in our study of the large-n asymptotics of C n (v) and Z n (q, v) at complex 
v [55]. 

3.4.2 Convolution formulae 

Specializing (3.40)-(3.43) and exploiting the symmetries of the complete graph, 
we obtain 

Z n (q 1 + q 2 ,v) = ^ (j^j Z k {q x , v) Z n - k (q 2 , v) (3.52) 

Z n (q 1 + q 2 ,v) = E ( k _ 1 J Z k(qi,v)Z n _ k (q 2 ,v) for n > 1 (3.53) 
k=l ^ ' 

Z n (q,v) = y~]l, 1 ) gCfc(^) Zn-k(q, v) forn>l (3.54) 
fc=i ^ ' 



A ._ 1 J(?i + ?2) - ^Jgi 



Z fc (gi, u) ^ n -fc(g2, u) for n > 1 

(3.55) 



On the other hand, these are just the one-dimensional identities (1.5), (1.6), (1.9) and 
(1.10) applied to the family a n (q) = Z n (q,v), which is indeed of the required form 
(1.2) by virtue of the generating-function formula (3.9). 
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The formula (3.54) can be used to compute the Z n (q,v) inductively given the 
Ck(v). On the other hand, if we specialize (3.54) to q — 1 and use (3.3), we obtain 
the identity 



(1 + v )»("-i)/2 
or equivalently 



it (fc _ |) ^ + ^) {n ~ fc)(n -^ 1)/2 for n > 1 (3.56) 



CJv) 



'1 + v 



,n(ra-l)/2 



-E(r 1 1 )^)(i+«) (B - fc,(B - 
fc=i ^ ^ 



■fc-l)/2 



for n > 1 , 
(3.57) 

which can be used to compute the C n (v) inductively ab initio. We call (3.56)/(3.57) 
the "linear" identity for C n (v). It goes back at least to Leroux [37, eq. (3.3)] and is 
probably much older; it can be proven by an easy direct argument based on consid- 
ering the size fc of the connected component containing a fixed vertex. [The same 
argument proves (3.54).] 

Alternatively, if we specialize (3.55) to qi — 1 and use (3.3), we obtain 



z n (q,v) = 



k=l 



n — 1 
fc- 1 



(1 + v )*(*-i)/2 Z n _ k (q,v) forra>l, 

(3.58) 

which allows us to go directly from from Z n (l, v) = (1 -fti)™*™ -1 )/ 2 to Z n (q, v) without 
passing through C n (v). We can also rewrite (3.58) as 



ra-l 



Z n (q,v) 



:i+u) n(n-l)/2 + J- 



k=l 



n — 1 
fc 



9 - 



n — 1 
fc - 1 



;i + u) (n - fc)(n - fc - 1)/2 z fc ( g ,v) 

for n > 1 , (3.59) 

which manifestly generalizes (3.57) and reduces to it when q — 0. 



3.5 A generalization 

Let a = {a n }^ be an arbitrary sequence of coefficients belonging to the ring R, 
satisfying ao — 1, and define the family of polynomials 

M 

Z n (g;a) = ]T q^\{a Ul \ (3.60) 

with the convention Z = 1. Note that Z n (l;a) = a n . This definition generalizes 
(3.47)/(3.51), and if we specialize to a n — (l + ■u) n < n ~ 1 )/ 2 we obtain Z n (q,v). 
A simple counting argument shows that 

n b / \ b 

*.<«»> = E£ E . . < 3 -«> 

b=0 ni,...,n 6 >l V 7 i=1 

E n i = n 
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Using this to compute the exponential generating function of the {Z n (q; a)}, we find, 
after a short calculation using the binomial series, that 



oo ,, 



nl 

n=0 




£-T^ a ) = £^«n • (3-62) 



In other words, the family {Z n (q;SL)} is of the form {a n (q)} defined in (1.2), with 
a n (l) = a n . But this should hardly be surprising, as (3.61) is simply the one- 
dimensional case of (2.12) with q 2 = q and q\ = 1. We have thus come full circle. 
Likewise, the recursion 



Z n (q;a) = q ^ K _ 1 J a k Z n - k (q - 1; a) 
k=i ^ ' 



(3.63) 



can be proven from (3.60) by an easy direct argument based on considering the size 
k of the block of 7r containing some fixed element of [n]. But this identity is nothing 
other than (1.11) specialized to q 1 — 1 and q 2 — q — 1. 



3.6 A related problem 

Suppose we try a definition analogous to (3.60), but with ordinary powers g' 71 "' in 
place of falling factorials q^. That is, let c = {c n }™ =1 be an arbitrary sequence of 
coefficients belonging to the ring R, and define the family of polynomials 

M 

Y n {q-c) = X> W n c W (3.64) 

7ren n i=l 

with the convention Yq — 1. Of course, the factor g' 71 "' is superfluous, because we can 
simply multiply each c n by q, but it is convenient to keep it explicit. This definition 
generalizes (3.46), and if we specialize to c n = C n (v) we obtain Z n (q,v). 
A simple counting argument shows that 

n ft , \ b 

*fa;c> = £ff e (,..::.,,) 11- (»■«) 

^— m, . . . ,nj, > 1 i 

Using this to compute the exponential generating function of the {Y n (q; c)}, we find, 
after a short calculation using the exponential series, that 

^ — F n (g;c) = expUJ] — c n . (3.66) 

n=0 ' \ n=l ' / 

In other words, the family {Y"„(g;c)} is again of the form {a n (q)} defined in (1.2). 
But this should hardly be surprising, as (3.65) is simply the one-dimensional case of 
(2.3) [i.e. (1.3).] We have again come full circle. 
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Likewise, the recursion 



n f n —\\ 

Y n (q;c) = g f J c fc y„_ fc (g; c) , (3.67) 
k=i ^ ' 

which can be proven from (3.64) by considering the size k of the block of n containing 
some fixed element of [n], is nothing other than (1.9). 

4 Nonlinear identity for the multivariate Tutte poly- 
nomial 

4.1 General identity 

In the preceding section we proved the identities (3.41) and (3.42) for the multi- 
variate Tutte polynomial Z G (q, v), which are based on choosing a single distinguished 
vertex i. We shall now prove a different identity that is based on choosing a pair of 
distinguished vertices Unlike the identities discussed in the preceding section, this 
one relies on the graphical structure of Z G (q, v) and does not appear to generalize to 
arbitrary families {a n (q)} of the type (1.2). 

So let G = (V, E) be a finite undirected graph with |V| > 2, and let i,jEV with 
% ^ j. We then have the following identity: 

Theorem 4.1 Let G = (V, E) be a finite graph, and let i,j e V with i ^ j. Then 



w cv 

W3i,Wjj 



q-l+ J] (l + O 

eeE(WJ) 



C G[w] (v) Z G[v \ w] (q, v) (4.1) 



where E(W,j) denotes the set of all edges with one endpoint in W and the other 
endpoint at j. In particular, 



Cq(v) = £ 



w cv 



n - 1 

eeE(W,j) 



C G [w](v)C G[v \ w] {v) . (4.2) 



We do not know whether there exists a generalization of (4.1) involving q\ and q2 
(rather than just q± = 0). 

The proof of Theorem 4.1 is based on looking at the connected component of i in 
the induced subgraph where j is deleted: 

PROOF. Start from the definition Z G (q,v) = q k{A) II v e- Let G' = (V, A) and 

ACE e£A 

G" = G'\j, and let W be the vertex set of the connected component of G" containing 
i. Let us now sum over all A C E that give rise in this way to a specified set W CV, 
and let us split this sum into two parts according as A does or does not contain at 
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least one edge from j to W. If A does not contain such an edge, then W is a connected 
component of G', giving rise to a factor q, and we get qCa\w\{y) Zg[v\w](q, v )- O n the 
other hand, if A contains at least one such edge, then W forms part of the connected 
component of j (e V\W) in G'; hence there is no factor q, but there is a factor v e for 
each edge e in the (nonempty) subset AnE(W,j). Summing over all such nonempty 
subsets of E(W,j), we get 



n ( i+v ^ - 1 

e&E(W,j) 



Cg[w](v) Za[v\w](q,v) . 



(4.3) 



Putting everything together gives (4.1). Specialization to q = yields (4.2). □ 

There is also a variant of (4.1) in which the right-hand side is completely expanded 
out. For notational simplicity we shall assume that G has no loop at the distinguished 
vertex j. (In the general case one must multiply by a factor 1 + v e for each loop e at 

./•) 

Proposition 4.2 Let G = (V, E) be a finite graph, let j G V , and suppose that G 
has no loop at j. Then 



Za(q,y) = E II 

n€U(V\{j}) Ben 



q-l+ (1 + 

eeE{B,j) 



Cg[b](v) 



(4.4) 



where the sum runs over (unordered) partitions n of V \ {j}, and the product runs 
over blocks B ofn. In particular, 



caw = e n 

7ren(v\{j}) Be-* 



n ( i+v ^ - 1 

eeE(B,j) 



Cg[b](v) 



(4.5) 



Proof. In (3.2), consider a term ACE, and let G' = (V, A) and G" = G'\j. Let 
7r be the partition of V \ {j} into vertex sets of connected components of G" . We 
can recover G' from G" by adjoining, for each block B of 7r, zero or more edges from 
the set E(B,j). If we adjoin zero edges, we get an extra factor q because B becomes 
the vertex set of a connected component of G' that is distinct from the component 
containing j; if we adjoin one or more edges, we get no such factor. This proves (4.4). 
Specialization to q = yields (4.5). □ 

When v e — v for all edges e, the formulae (4.4) and (4.5) can be found in Gessel [17, 
Theorems 13 and 10]. The special case of (4.5) in which v e = —1 for all e was also 
found by Borgs [8] and used recently by him [9, Lemma 3.2] to bound the complex 
zeros of chromatic polynomials (a variant of the proof in [56]). Indeed, we were led to 
formulate (4.4)/(4.5) by meditating on Borgs' special case, oblivious to the fact that 
they had already been essentially found by Gessel! 
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Finally, let us use (4.2) to prove an interesting inequality concerning the polyno- 
mials C G (v). Let 

c(A) = \A\ - \V\ + k(A) (4.6) 

be the cyclomatic number of the subgraph (V, A), and let us define the generalized 
connected sum 

C G (v,A) = J2 XC(A) Hv e (4.7a) 

ACE eeA 
k{A) = 1 

= A- (|l/| - 1} C G (Av) . (4.7b) 

Of course, (4.7b) shows that C G (v, A) contains no more information than C G (v); it is 
simply a convenient way of scaling all the variables v e simultaneously while removing 
a factor A' v ' _1 . In particular, C G (v, A) interpolates between the spanning-tree sum 
(A = 0) and the connected-spanning-subgraph sum (A = 1). We then have the 
following result [56, Remark 2 in Section 4.1] [53, Proposition 2.5]: 

Proposition 4.3 Let G = (V, E) be a finite graph (V ^ 0) equipped with real edge 
weights v = {y e ) e& E satisfying — 1 < v e < for all e G E. Then 

(-l)' +|yhl J^ G (v,A) > (4.8) 
on < A < 1, for all integers t > 0. 

This inequality was proven in [53,56] using a "partitionability" method going back 
to Penrose [41]. 17 Here we prove it using the recursion (4.2): 

Proof of Proposition 4.3. By induction on \V\. When \V\ = 1, G necessarily 
consists of zero or more loops attached to the sole vertex, so 

C G (v,A) = H(l + Xv e ). (4.9) 

Using v e < and 1 + Xv e > 0, it is easy to see that (4.8) holds. 

Now assume that \V\ > 2. Replacing v by Av in (4.2) and using (4.7b), we obtain 



C G (v,A) = Yl A_1 



w cv 

W Bi,WJ j 



n - 1 

e€E(W,j) 



CG[W](y, A) C G [ V \ W ](v, A) 



(4.10a) 



E E ^ Bhl U v e) c G[w](y,x)c G[ v\w](y,x)- 

W CV \<2^BCE(W,j) e€B 
W 3i,W $ j 



(4.10b) 



17 For the partitionability method (which applies to matroids as well as graphs), see also [7,19] 
and the other references mentioned in [53, Section 2.2]. 
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Now apply cr/d\ to both sides: on the right-hand side we will have terms in which 
i\ derivatives act on A' B ' _1 Y[ v e, ^2 act on Cc[w] ( v j A) , and £ 3 act on Cc[v\w] ( v j A) , 
where ^1 + ^2 + ^3 = ^- The induction hypothesis is applicable because 1 < \W\ < 
|Vj — 1. So it suffices to check that 

("D'-'^rf £ ^'"Th) so. (4.11) 

To prove this, let us order the elements of E(W,j), i.e. consider E(W, j) ~ [m] where 
m = \E(W,j)\. Then, by singling out the smallest element in each set B, we can 
write 

m m 

e ^ M n^ = i> n^+M). (4.12) 

0^BC[m] e€B i=l j=i+l 

We have Vi < 0; each d/d\ brings down some factor Vj < 0; and the undifferentiated 
factors satisfy 1 + Xvj > 0. So (4.11) holds. □ 

It is worth remarking that (4.11) is nothing other than Proposition 4.3 specialized 
to the graph consisting of two vertices connected by m parallel edges. Moreover, 
the identity (4.12) is nothing other than (one version of) the partitionability identity 
for Ki m) . 

4.2 Specialization to the complete graphs 

Specializing (4.1) to the complete graph K n with equal weights v, we obtain 

n 1 / — 2\ 

Z n (q, v) = [ k _ ! J [9 + (! + v ) k ~ !] C ^ v ) v) forn > 2 (4.13) 

(see [17, eq. (5)] for an equivalent formula). Specializing this to q — 0, we obtain 

Cn(v) = ( I _ 1 ) K 1 + v ) k ~ !] C'*^) C -k( v ) forn > 2 , (4.14) 
fc=i ^ ' 

a result apparently first proven by Leroux [37, eq. (3.5)] (see also [17, eq. (2)]). 18 And 
specializing this latter formula to v — 1, we obtain the identity of Riordan, Nijenhuis, 
Wilf and Kreweras [31,40,47] for counting connected graphs. We call (4.14) the 
"nonlinear" identity for C n (v). 

On the other hand, specializing (4.13) to q — 1, dividing both sides by (1 + v ) n ~ 1 
and relabelling n — > n + 1, we recover the "linear" identity (3.56). 

The formula (4.14) has also been derived in a different context. Let T be a tree 
with vertex set [n], rooted at the vertex 1. An inversion of T is an ordered pair (j, k) 

18 See also [5, p. 306, Exercise 4.2.2] for a proof of (4.14) in species language. 
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of vertices such that j > k > 1 and k is a descendant of j (i.e., the path from 1 to 
k passes through j). We define the inversion enumerator for trees [20,31,38] (see 
also [10, 17, 19]) to be the polynomial 

Uv) = E y im(T) ( 4 - 15 ) 

trees T on [n] 

where inv(T) denotes the number of inversions in T. This polynomial turns out to 
be related to C n (v) by the beautiful formula [3,20,31,38] 

C n (v) = v n - l I n (l+v). (4.16) 

Now, Mallows, Riordan and Kreweras [31,38] show that I n {y) satisfies the recursion 

J M = E (l Z f) (XV ) J M 7 »-*(v) forn > 2 . (4.17) 
k=i ^ ' \j=o J 

But using (4.16), it is easily seen that (4.17) is equivalent to (4.14). 

Let us remark, finally, that even without using (4.15), it follows immediately 
from the recursion (4.17) [together with the initial condition h(y) = 1] that the 
polynomials I n (y) = C n (y — l)/{y — l) ra ~ 1 have nonnegative (indeed, strictly positive) 
integer coefficients. On the other hand, the nonnegativity of the derivatives of I n (y) 
at y = is also a special case of Proposition 4.3. 



A A generalization of Mobius inversion on the par- 
tition lattice 

In this appendix we present an apparently new generalization of Mobius inver- 
sion on the lattice of partitions of a finite set, which is inspired by the formulae in 
Section 3.3.1 and in particular by (3.30), and more generally by (2.13). We refer 
to [59, Chapter 3] for basic facts about posets and Mobius inversion. 

To begin with, let P be a finite poset. Then the zeta function ( on P is the 
function ( : P x P — > Z defined by 



C(x,y) 



1 if x < y 
if x ^ y 



(A.l) 



(We think of ( as a matrix whose rows and columns are indexed by P.) The Mobius 
function /i on P is the two-sided matrix inverse of (, i.e. it satisfies (fi = fi( = I. It 
can be computed by the recursion 



/i{x,y) 



1 



if x 



y 



- E »(x,z) iix<y 

x<z<y 

if x % y 



(A.2) 
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Now let S be a finite set, and let P be the lattice H(S) of partitions of S, ordered 
by refinement. It is well known [59, p. 128] that the Mobius function of H(S) is given 
by 



y](_l)K-i(\. _ i)i if ° - n = i B i> ■■■^ B k] and Bi is 
, M partitioned into A, blocks in a . . „. 

M^) = < ( A -3) 

if (7 ^ 7T 

We generalize this as follows: let gi and q 2 be indeterminates, and define 

{A Yr^ga - ?gi) if a ~ 71 = ■■■' Bk ^ and Bl is 

/Jj /-A partitioned into A« blocks in a . . , . 

.-i j-i ^ ( A .4) 

if (7 £ 7T 

Let us observe for future reference that 

m— 1 

= ^ 2 _1 feM)- = (zr'Wft - !)— • ( A -5) 

J'=l 

Then we have in particular 

A*o,i = C (A.6a) 

A*i,o = (A.6b) 

IJ, qjq = I for all q (A. 6c) 

We shall prove the following generalization of Mobius inversion: 

Theorem A.l Let S be a finite set, and define the matrices [i q ^ m on U(S) by (A.4). 
Then 

/V,«2 ^qim = /V,?3 (matrix multiplication) (A. 7) 

We interpret (A. 7) as an identity in the polynomial ring q 2 , q^], but of course it 
also holds when q±,q 2 , ?3 are specialized to specific integer values (or more generally to 
specific values in a commutative ring R). In particular, Mobius inversion corresponds 
to the special cases (q 1: q 2 , (fo) = (0,1, 0) and (1, 0, 1). 

The proof of Theorem A.l will be based on the following lemma, which we think 
is of some interest in its own right. It is probably not new, but we have been unable 
to find any reference. 

Lemma A. 2 Let m be a positive integer, and let r and s be indeterminates. Then 

J2 rM| JJ sM ) = M- ( A - 8 ) 
as an identity in the polynomial ring Z[r, s] . 
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Let us remark that if one divides (A. 8) by s m and takes s — > oo (or equivalently 
just extracts the coefficient of s m ), one obtains the well-known formula (3.21)/(3.38); 
while if one sets r = q/e, s = eq and takes e — > 0, then one obtains the formula 

m(6,tt) gl-l = ^ (A.9) 

7ren m 

due to Rota [50, section 9] [59, pp. 128, 162, 187] (here denotes the partition in 
which every element is a singleton). 

Proof of Lemma A. 2. It suffices to prove (A. 8) for positive integers r, s. Then 
the right-hand side counts the proper colorings of the complete graph K m with rs 
colors. But so does the left-hand side, if we use the color set [r] x [s] and define u to 
be the partition of [m] in which two vertices are placed in the same block if and only 
if they receive a color with the same first index. □ 



Second Proof of Lemma A. 2 (suggested independently by Christian Kratten- 
thaler and Richard Stanley). We use the isomorphism of weighted species 

Partitions = Sets o NonemptySets , (A. 10) 

where a Set of cardinality b is given a weight r-, a NonemptySet of cardinality k is 
given a weight s-, and a partition uj is weighted as on the left-hand side of (A. 8). 19 
The corresponding exponential generating functions are 

F(x) = E^r* = (l + *) r ( A - n ) 

6=0 

OO 

G ^ = Ex,*" = (! + *)*-! (A.12) 

k=l 



and hence 



oo _ 
X r 



(FoG)(x) = (l+*r = E-fM-' (A.13) 
which proves (A. 8) for all m > 0. □ 



mi 

m=0 



Proof of Theorem A.l. Let a < n = {B ± , . . . , B k } where B { is partitioned into 
Aj blocks in a, and let us compute 

(/W ^92,93) T) = E Vquq2(<7,T)Hq 2iq3 (T,ir) . (A.14) 

(7<T<TT 



19 See [5, pp. 44, 86] for this isomorphism with a different weighting. For the reader unfamiliar with 
the theory of combinatorial species [5] , it suffices to observe that a partition of a finite set S is simply 
a set of nonempty subsets of S that disjointly cover 5; one can then invoke [60, Theorem 5.1.4] to 
complete the proof. 
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The partition r is specified by saying, for each % e [&], how the corresponding \ blocks 
of a get grouped in r. The sum (A. 14) will then factorize over i, with the ith factor 
given by 

(Mi |B a |-l \ in-l 

n n fe-^i) n^-M 
a=l j=l J k=l 

u = {B 1 ,..:,B IH } 

E ( I! 9P a 'fe '(ft/ft)^ 1 ) <#<fe '(fc/fc)* 

w = {Bi,...!b w } 

= siV 1 E ( ftteM)^) 

u> e n A . \a=i / 

u = {B 1 ,... t ,B IH } 

= qfq^ite/qi)^ (A.15) 

by Lemma A. 2. □ 



Second Proof of Theorem A.l (suggested by Richard Stanley). We prove (A. 7) 
for S = [n] simultaneously for all n > 1, using the interpretation of convolution of 
multiplicative functions on II = (n 1 ,Il2, . . .) in terms of composition of exponential 
generating functions [60, pp. 7-8]. We first observe that n qxm is a multiplicative 
function on II, with underlying numerical function 

m— 1 

/»•*("») = II = ^fcW?!)" (A.16) 

j'=i 

and corresponding exponential generating function 

oo m 

W*) = E /<zi.«( m ) = ?2 _1 [(1 + ?i^) 92/91 " 1] • (A.17) 

Til. 

m=l 

It is easy to check that 

(F quq2 (x)) = F qim (x) . (A.18) 

By [60, Theorem 5.1.11], the identities (A. 7) for all n > 1 are an immediate conse- 
quence. 20 □ 



20 Using [60, Theorem 5.1.11] with the coefficient field K taken to be (say) K., this argument 
proves (A. 7) for 91 , (fe , <?3 £ K \ {0}, which is more than enough to imply it as a polynomial identity. 
Alternatively, we can use [60, Theorem 5.1.11] with K taken to be (say) the field M.(qi, q2, (73) of 
rational functions in the indctcrminatcs qi,q2, <?3, which yields (A. 7) directly as a polynomial identity. 
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Remark. If we define the diagonal matrix D r by 

f r H if cr = 7T 
DrM = \ (A.19) 
10 it a f: ix 

then an easy calculation shows that 

D r fJ, qi i?2 D r = flrq 1 ,rq 2 ■ (A. 20) 

This is why many things involving fJ, quq2 depend only on the ratio 92/91 [cf. (3.28)]. 
□ 



We can apply Theorem A.l as follows: Let R be a commutative ring, and fix 
a function /: H(S) — > R and an element go £ R (it could be an indeterminate if 
desired). Now introduce a new indeterminate q, and define F q : H(S) — > by 

Fq = f^q ,q (A.21) 

or in more detail 

W = E /(^)/V>^)- (A.22) 

aGn(5) 

By construction we have imposed the "initial condition" 

F qo = f (A.23) 

(here F qo means of course that the indeterminate q is replaced by the value qo E R). 
Most importantly, Theorem A.l implies that the (/-dependence of F q is "coherent" in 
the sense that 

F qi H qi , q2 — F q2 . (A. 24) 

The formula (3.30) is now the special case of (A. 24) in which we take S — V, 

F q (7r) = J] Z G[B] (q,v) , (A.25) 

and we evaluate (A. 24) at n = 1 (the partition in which all of V belongs to a single 
block). More generally, in place of Zg[b](q,v) we can use <2i s (<?) where {a n (g)} is any 
family arising as in Section 2; we then obtain (2.13) specialized to n = l. 21 

The formulae in this appendix seem to be special for the partition lattice. Still, it 
may not be totally absurd to ask: Might there exist u q\-q2 generalizations" of Mobius 
inversion for some other posets? Maybe even for some natural class of posets? 
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Note also that, in these formulae, we use a n (q) [or c n ] only for < n < 1. 
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